A c / A c production asymmetries in pp and tt p 

collisions 

G. Herrera * 

Centro de Investigacion y de Estudios Avanzados 
Apdo. Postal 14-740, Mexico 07000, DF, Mexico 

J. Magnin ' 

Centro Brasileiro de Pesquisas Ffsicas 
Rua Dr. Xavier Sigaud 150, CEP 22290-180, Rio de Janeiro, Brazil 



PACS 13.60.Rj (Baryon production) 

13.87.Fh (Fragmentation into hadrons) 
14.20.Lq (Charmed baryons) 

Abstract 

We study A c /A c production asymmetries in pp and ir~p collisions 
using a recently proposed two component model. The model includes 
heavy baryon production by the usual mechanism of parton fusion and 
fragmentation plus recombination of valence and sea quarks from the 
beam and target hadrons. We compare our results with experimental 
data on asymmetries measured recently. 
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1 Introduction 



Recent experimental data on charmed meson hadroproduction shows a strong 
correlation between the flavor content of the incident hadron and the pro- 
duced meson |], 0. This effect, known as leading particle effect, has been 
observed also in A c production in pp [|, f| and n~p || interactions, in A& pro- 
duction in pp collisions and S c production in hyperon-nucleus interactions 

1- 

The leading particle effect can not be explained within the usual frame- 
work based in the factorization theorem. In fact, although perturbative QCD 
(pQCD) at Next to Leading Order (NLO) predicts a small enhancement of 
the c over the c-quark cross section || , the effect of this asymmetry is very 
small to account for the observed D / D + asymmetry measured in 7r~-nucleus 
interactions 0. On the other hand, the evidence of leading particle effects in 
A c production at large xf clearly shows that non-perturbative mechanisms 
play a fundamental role in hadroproduction. 

Leading particle effects in charmed hadron production has been studied 
from two different points of view: the soft scattering - low virtuality approach 
JJ, in which charmed meson anti- meson production is reviewed within the 
valon recombination model, and two component models, in which charmed 
hadron production takes place by two different mechanisms, namely pertur- 
bative QCD (pQCD) followed by fragmentation plus coalescence of intrinsic 
charm [[H]] or recombination of the valence quarks with the charmed sea 



quarks of the initial hadrons liberated in the collision [fLl [. In the two com- 
ponent models, the intrinsic charm coalescence or the valence and sea quark 
recombination gives the non-perturbative contribution which takes into ac- 
count the observed flavor correlation. 

Leading particle correlations can be quantified by studies of the produc- 
tion asymmetries between leading and non-leading particles. The production 
asymmetry is defined by 

A ( x ) — d(J L/dx F - da NL /dx F 
L/WL F d<jL./dxF + daNi/dxF 

in which L stands for Leading and NL for Non-Leading particles. In this 
paper we give the predictions of the recombination two component model for 
A c /A c production asymmetries in pp and ir~p interactions. We also compare 
the predictions of the model with recent experimental data in pp 0] and n~p 
[HI interactions. 
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2 A c and A c production in pp and n p collisions 



In this section we review A c /A c production in the recombination two com- 
ponent model [11 1 . In section 



J] we present the calculation of the parton 
fusion contribution and in section B72] the recombination picture is studied 
for both pp and ir~p interactions. 



2.1 A c and A c production via parton fusion 

The calculation we present here is at Lowest Order (LO) in a s . A constant 
factor K ~ 2 — 3 is included in the parton fusion cross section to take into 
account Next to Leading Order (NLO) contributions |I2[ . 

The inclusive xp(— 2pi/y/s) distribution for A c (A c ) production in hadron- 
hadron interactions, assuming factorization, has the form [ |T3|j 



da^ 

dXF 



^Vs J H^ b B (x a ,x b ,Q 2 )^ DAc/c Z dzdp 2 T dy , 



(2) 



where 



H^ b B (x a ,x b ,Q 2 ) = T Jafi [q a (x a ,Q 2 )q b (x b ,Q' 2 



+ q a (x a , Q 2 )qb(x b , Q 2 )) S \ qq - 

' dt 

+9a(x a ,Q 2 )g b (x b ,Q 2 )- I \ g g , 

dt 



(3) 



with x a and x b the parton momentum fractions in the inital hadrons A and B, 
q(x, Q 2 ) and g(x, Q 2 ) the quark and gluon distribution in the corresponding 
colliding hadrons, E the energy of the produced c (c)-quark and -Da c / c (z) the 
appropriated fragmentation function. In eq. ||, p\ is the squared transverse 
momentum of the produced c (c)-quark, y is the rapidity of the c (c) quark 
and z = xpjx c is the momentum fraction of the charm quark carried by the 
A C (A C ). The sum in eq. ||| runs over a,b = u,u,d,d,s,s. da/di \ q q and 

da/dt \ gg are the elementary cross sections for the hard processes qq — > cc 
and gg — > cc at LO given by 



da 
dt 



ira; 



(Q 2 ) cosh (Ay) + m 2 /rh. 



i<?<?~ 



9m^ 



da 
dt 



ira; 



[1 + cosh (Ay)Y 



(Q 2 ) 8cosh(Ay)-l 



199" 



96m4 



[1 + cosh (Ay)] ,: 



* / a \ 2m 2 
cosh (Ay) + 



mi 



+ 



1m\ 
mi 



(4) 
(5) 
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where Ay is the rapidity gap between the produced c and c quarks and 
rh 2 = m 2 + p\. The Feynman diagrams involved in the calculation of eqs. ^ 
and |5| are shown in Fig. [TJ 

For consistency with the LO calculation, we use the GRV-LO parton dis- 



tributions for both the proton [|TJ] and the pion [15|. The hard momentum 
scale is fixed at Q 2 = 2m 2 . 

The fragmentation is modeled by two different functions; the Peterson frag- 



mentation function extracted from data in e + e interactions 16 



N 

D A c /c{z) = — — T72 (6) 

Z[l-1/Z- E c /(l- Z)\ 

with e c = 0.06 and the normalization defined by J2h I Du/c{ z )dz = 1 and 
the delta fragmentation function 

D Ac/c (z) = 5(l-z). (7) 

The use of the delta fragmentation function implies that the A c is produced 
with the same momentum carried by the fragmenting c-quark. This mech- 
anism for fragmentation has been used to simulate the coalescence of the 
c-quark, produced in a hard interaction, with light valence quarks coming 



from the initial hadrons ITU 



In Figs. and ^ we show the parton fusion distributions obtained for the 
two fragmentation functions in pp and ir~p interactions respectively. Notice 
that the A c and the A c parton fusion distributions are equal at LO. However, 
at NLO there is an small c/c asymmetry due to the fact that to this order 
in perturbative theory the cross section for the production of a quark differs 
from the cross section for the production of an anti-quark || and consequently 
a small A c /A c asymmetry is present. 

2.2 A c and A c production via recombination in pp and 
7r~p interactions 

The method introduced by K.P. Das and R.C. Hwa for meson recombination 



was extended by J. Ranft |Tj| to describe single particle distributions of 
baryons in pp collisions. In a recent work R.C. Hwa || studied -D 1 * 1 inclusive 
production using the valon recombination model. 

In recombination, all products of the reaction appearing in the forward 
region (xp > 0) are thought as coming from the beam particle fragmentation 
while hadrons produced in the backward region {xp < 0) come from the 
target fragmentation. Since the two reactions we are considering here have 
the same target particles, then the recombination A c and A c inclusive xp 
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distributions are the same for both reactions in the backward region. In the 
forward region, however, they must be different since the outgoing hadron is 
produced from the debris of different initial particles. Hence it is sufficient 
to study A c and A c production by recombination in the forward region for 
pp and ir~p collisions, being the backward Xp distributions symmetrical with 
respect to the forward xp distribution in pp interacctions in both reactions 
under study. 

The invariant xf inclusive distribution for A c (A c ) production in pp 
interactions in the forward region is given by 

2E da rec [ x f dxi dx 2 dx 3 a c (a c ) , 
rer 7= -; =/ ^3 {xi,X2,x 3 )R 3 {x 1 ,x 2 ,x 3 ,x F ) (8) 

<7 rec ^S dXF Jo X\ X2 X 3 

where E is the energy of the produced hadron, a/s is the center of mass en- 
ergy, R 3 (xi,X2, x 3 , xf) is the recombination function and F 3 ° ^ c ^ (xx,x 2 , x 3 ) 
is the three quark distribution function. Xi, i=l, 2, 3 is the momentum frac- 
tion of the i th quark with respect to the proton. 

Following the approach of Ref . [ , the three quark distribution function 
is assumed to be of the form 



F 3 Ac(Ac) (xi, x 2 , x 3 ) = (3 p Fl ( Xl ) Fl (x 2 ) F£ (x 3 ) (1 - x x - x 2 - x 3 f* (9) 

In eq. |||, F Qi = Xiqi(xi) is the single quark distribution of the i th valence 
quark in the produced A c (A c ) inside the proton and the coefficients (3 P and 
7 p are fixed using the consistency condition 

F q (xi) = / dxj / dx k F 3 c{ c > (x x ,x 2 ,x 3 ) . 

Jo Jo 

i, j, k = 1,2, 3 (10) 
which must be valid for the valence quarks in the proton. 



In eq. || we use the single quark GRV-LO distributions [|Tj]. Since the 
GRV-LO distributions are functions of the momentum fraction x and the 
momentum scale Q 2 , then our F^ c ^ Kc {x\,x 2 , x 3 ) also depends on Q 2 . We use 
Q 2 = Am 2 c with m c = 1.5 GeV as in Ref. |IT]. Notice that the scale in the 
recombination is fixed at a different value than in parton fusion. In fact, in 
parton fusion the scale Q 2 is fixed at the vertices of the Feynman diagrams 
involved in the perturbative part of the calculation while in recombination 
this parameter must be chosen in such a way that small variations in its value 
does not change appreciably the charmed sea of the initial hadron. 

Concerning the single-quark contributions to the three quark distribution 
function, some comments are in order. In pp interactions, contributions to 



5 



the A c inclusive xp distribution in the large xf region come mainly from 
Valence- Valence-Sea (VVS) recombination processes, being other processes 
involving more than one sea-flavor in the recombination completely negligible 
due to the fast fall of the sea-quark distributions. Conversely, in the small 
xp region, contributions of Valence-Sea-Sea (VSS) and Sea-Sea-Sea (SSS) re- 
combination processes are important since sea-quark distributions are peaked 
at xp about zero. In this region, the later processes dominate largely over 
VVS recombination. Hence, in order to have a more acurate prediction for 
the A c /A c asymetry at small xp, VSS and SSS processes must be included 
in the calculation of the A c cross section. A c production in pp interactions 
proceeds only through SSS recombination. 
For the recombination function we take 

(XiX2) ni X n2 

R 3 {xx,x 2 , x 3 ) = a — ni+n2 _i $ (xi +x 2 + x 3 - x F ) (11) 

Xp 

allowing in this way for a different weight for the heavy c (c) quark (marked 
with index 3) than for the light u (u) and d (d) quarks (indexed 1 and 
2 respectively), as suggested in Ref. in conection with charmed meson 
production in the valon model. 



The constant a in eq. |TT| is fixed by the condition [11 



i da rec 
dx F - = 1 (12) 

dXF 

for A c inclusive production. Using the same value for the parameter a ob- 
tained in eq. [12] to normalize the A c xf distribution, the anti-particle cross 
section is given relative to the particle cross section. In this way a r ^ has 
the physical interpretation of the A c total recombination cross section and it 
should be fixed by experimental data. 

In Fig. [|we show the A c and A c xp inclusive distribution for n\ — n 2 = 1 



lq] and n\ = 1, n 2 — 5 Q with f3 p = 75 and 7 P = —0.1 as given in Ref. [[11 

In 7i~p interactions in the beam fragmentation region both A c and A c 
inclusive xp distributions are given by formulas formally identical to that 
of eq. H As the 7r~ has an u and a d valence quarks, the recombination 
processes involved in A c as well as in A c production are VSS and SSS. On 



the other hand, since the GRV-LO [15 valence distributions and quark and 



anti-quark sea distributions in the pion are equal, then the A c and A c inclusive 
xf distributions are the same. Notice that A c production in ir~p interactions 
in the proton fragmentation region is double leading and A c is non-leading 
whereas in the beam fragmentation region both particle and anti-particle 
production are leading, so in this A c /A c asymmetry is expected for 

xp < 0. 
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3 A c /A c production asymmetry in the recom- 
bination two component model 

The total inclusive x F distributions are obtained by adding the parton fusion 
and recombination contributions given by eqs. [| and § respectively. Its 
general form is 



da tot da pf rec da rec 

+ cr, 



dxi 



W " dx F ^" Ac dx f 



Ac(Ac) 



(13) 



where we have uncovered the parameter a r ^ hidden in the da rec jdxp def- 
inition of eq. H Let us stress that o r ^ is the only free parameter in the 
model. 

The A c /A c production asymmetry is obtained by replacing the corre- 
sponding total inclusive distribution given by eq. [13] in eq. [I]. In this way 
we obtain 

A , \ rec da rec /dx F \ Ac -da rec /dx F | x 

/a \ x f) = cr A f t (14) 

Ac/Ac c 2daPf/dx F + a r A e c c \ • 

for the asymmetry function of x F . 

It should be noted that pQCD at NLO predicts a bigger A c than A 
production at large values of x F . Since the asymmetry coming from pQCD 
is of the order of 10 % at x F ~ 1 and in the large x F region the recombination 
component is several orders of magnitude bigger than the perturbative part, 
this contribution to the asymmetry is completely negligible. 

In 7i~p interactions in the forward region (x F > 0), the A c and the A c re- 
combination cross sections are equal, then our model predicts zero or a small 
negative asymmetry due to NLO pQCD effects. This prediction is consistent 
with the ratio N(A C )/N(K C ) = 0.99 ± 0.16 measured by the ACCMOR col- 
laboration [Tj| which indicates no asymmetry (A = —0.002 ± 0.16). Also a 



c 



preliminary analysis in 500 GeV/c 7r - -nucleus interactions from the E791 col- 
laboration H is consistent with zero asymmetry in the region < x F < 0.55. 

In the backward region, the same E791 preliminary analysis shows a large 
asymmetry, reaching a value of A w 0.4 at x F ~ —0.13 but with large error 
bars §. The ratio A^(A C )/A^(A C ) = 1.17±0.08 (A = 0.078±0.034) measured 
by the E791 collaboration || is described by our model with a recombination 
cross section a A e c c ~ l.lcr p ^ and n\ — 1, n-i — 5 in the recombination function 
ofeq. |Tl"| in the Peterson fragmentation scheme. Here o v $ = Jq dx F da p ^ ' j 'dxF- 



Both Peterson and Delta fragmentation in combination with n% = ri2 = 1 in 
the recombination function tend to give a more slowly growing asymmetry 
with x F . 
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The E769 experiment H] found a lower limit of A = 0.6 in the forward 
region with a 250 GeV/c incident ir~ , and p beam on a multifold target of 
Be, Cu, Al and W. However, altough this experiment presents an evidence 
of A c / A c asymmetry, the value quoted in Ref . [|J for the asymmetry is not 
significative due to its very low statistic. 

In Figs. [5] and |6] we show the asymmetry as a function of xp predicted by 
the recombination two component model in pp and K~p interactions respec- 
tively with the recombination normalization suggested by the preliminary 
analisys from the E791 collaboration. 

4 Conclusions 

In this work we have presented the predictions of the recombination two 
component model for the A c /A c production asymmetry in pp and ir~p inter- 
actions. 

The model seems to describe adequately the scarce data available on 
the subject but, of course, more experimental data are nedeed not only on 
A c /A c asymmetries as well as on A c and A c xp distributions in order to do a 
meaningful comparison. 

Using the preliminary analysis of the E791 collaboration ||, the model 
predicts that the recombination cross section is as big as the parton fusion 
cross section, giving a clear idea about the extent to which non-perturbative 
contributions are important in charm hadroproduction. 

It should be noted that if the production asymmetry is independent of 
the energy, as indicated by measurements of charmed meson/ anti-meson 
asymmetries, the model predicts that the ratio of the recombination to the 
parton fusion cross section does not depend on the energy and could bring 
about important information about the non-perturbative processes in heavy 
hadron production. 

In the recombination two component model, the part responsible for the 
asymmetry is recombination. It involves different kinds of processes for A c 
than for A c production, giving xp inclusive distributions which are different 
in shape for the particle than for the anti-particle. 

Notice that strong diquark effects are present in the model for A c pro- 
duction in the proton fragmentation region in the two reactions considered 
in this paper. In fact, since the u and d are valence quarks in the proton, 
the ud diquark carries a large amount of the proton momentum when it is 
liberated in the collision. This large momentum is then transferred to the 
A c when the ud diquark recombines with a sea c-quark to form the outgoing 
particle. 
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Another model which has been used to make predictions on the A c /A c 
production asymmetry is the intrinsic charm two component model flQfl . In 
this model the enhancement in the A c over the A c cross section is due to the 
coalescence of u, d and c quarks coming from the \uudcc) Fock state of the 
proton. 

However there are important differences between predictions obtained 
with one or another model. 

In fact, one of the possibly most interesting features which distinguish the 
recombination from the intrinsic charm two component model is the ability 
of the first to produce large asymetries in the small xf region. This diference 
between the two models is due to the fact that recombination of sea and 
valence quark produces a A c xp distribution peaked close to zero whereas 
the intrinsic charm model gives a A c ^ distribution slowly growing from 
zero at xf = to its maximum at xf about 0.6. 

Finally, it should be noted that the two models seem to describe ade- 
quately the shape of the A c inclusive xf distribution (see Refs. [ 10|, [TTH ) 



but both models predict very different forms for the A c /A c asymmetry, so 
a meaningful comparison between models and experimental data should be 
done with data on asymmetries as well as on if distributions. 
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Figure Captions 



Fig. 1: Feynman diagrams involved in the LO calculation of the parton fusion 
cross section. 

Fig. 2: Parton fusion cross section for pp at 500 GeV/c energy beam. Peterson 
(dashed line) and Delta (full line) fragmentation are shown. 

Fig. 3: Parton fusion cross section for n~p at 500 GeV/c energy beam. Peter- 
son (dashed line) and Delta (full line) fragmentation are shown. 

Fig. 4: Recombination cross section in pp collisions for ri\ — 1, n 2 = 5 (full 
line) and m — n 2 — 1 (dashed line). Curves marked "SSS" are the 
corresponding A c distributions. The A c distribution is normalized to 
unity and the A c distribution is normalized relative to the first one. 

Fig. 5: A c /A c asymmetries in pp interactions for Peterson fragmentation and 
rii = 1, n 2 = 5 in the recombination function calculated for the recom- 
bination cross section suggested by the E791 preliminary data. The 
lower limit obtained by the E769 (Ref. [4]) experiment is shown. 

Fig. 6: Same as in Fig. 1 for n~p interactions. Experimental points are calcu- 
lated from the N(A C )/N(A C ) ratios measured in Refs. [5] and [19] 
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Figure 2: 
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Figure 4: 
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